Quasiparticle tunneling: P(E) theory 
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A calculation of the energy decay rate of a Josephson qubit from non-equilibrium quasiparticles 
is made using the environmental P(E) theory. For a large-capacitance qubit, we extend the theory 
to include the tunneling of quasiparticles, which has an electron- and hole-like charge components. 

PACS numbers: 



In this note, we calculate the dissipation due to the 
tunneling of quasiparticles using the environmental P(E) 
theory [1] to properly account for coupling to the qubit 
from the charge transfer. We will first calculate the 
damping rate for a inductor-capacitor resonator using 
environmental theory; the prediction for phase qubits, 
for example, will be similar since qubit and harmonic 
oscillator matrix elements are nearly identical. We will 
then consider the tunneling of quasiparticles, which has 
tunneling of charge carriers that is both electron- and 
hole-like. 

Standard environmental theory gives a probability 
P{E) for the environment to absorb energy from a tun- 
neling event with charge transfer q 
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where Z(u>) is the environmental impedance, Rk = h/q 2 , 
and we have assumed for simplicity T = 0. We first 
consider an inductor in parallel with a Josephson junc- 
tion of capacitance C, which has a resonance frequency 
tu r = 1/ y/LC. Here L is the parallel combination of the 
external inductance and the Josephson inductance. The 
impedance of the environment is 

Z(u) = -^[8{u + L0 r ) + 5{w - u r )\ , (3) 

which describes how the resonator can both absorb or 
emit photons at frequency ui r . Defining the resonator 
impedance Z r = l/u r C and inserting Eq. (3) into Eq. (2), 
we find 
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For the case of large capacitance resonators where Z r <C 
Rk, the second exponential can be expanded in a Tay- 
lor's series giving 



P(E) ~ (1 - irZ r /R K )6{E) + (irZ r /R K )5(E - tvaj r ) 



The probability for the environment to absorb a photon 
is given by the second term, and can be rewritten as 



p(hu r ) = 
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This probability is much smaller than one, which makes 
physical sense since the charging energy q 2 /2C of the 
tunneling event is much smaller than the environmen- 
tal energy hw r . For the case where the environment has 
initially one photon, the probability p{— oj r ) for the envi- 
ronment to emit a photon during tunneling has the same 
magnitude as Eq. (6). 

We have assumed in this calculation that charge q = 
—e is transferred across the junction in a tunneling event. 
This assumption must be modified since quasiparticles 
are both electron (q = — e) and hole (q = e) like, and their 
interference introduces coherence factors to the total tun- 
neling rate. The environmental theory must be derived 
in a way that properly accounts both for the quasiparti- 
cle states and the transfer of charge to the environment. 
The tunneling Hamiltonian is given by 
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where t LR is the tunneling matrix element. Changes 
to the internal (uncharged) state of the superconductor 
is described by the creation and annihilation operators. 
The L and R indices describe states in the left and right 
superconductor, and a describes the two spin states of 
the electron. The change in the electrical circuit to a 
tunneling event is represented by the charge displacement 
operators e ±v> , which describe the transfer of ±e from the 
tunneling event. Here, tp is the conjugate coordinate to 
charge q and corresponds to a dimensionless flux with 
dip/dt — (e/h)V = (e/h)(q/C). It has a commutation 
relation [q, ip] = ie and gives a charge displacement oper- 
ator according to the relation e lip qe~ l{p = q — e. 

The electron operators need to be re-expressed in terms 
of the quasiparticle operators 7 because quasiparticles 
are eigenstates of the superconducting leads. The four 
electron operators are given by 

Cfc| = Ufe7fcot + Vkltii Cfe + = u klkai - Vklti^ 
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The standard electron and hole occupation factors u k and 
Vk are given by 



u| = (l/2)(l + 
i£ = (l/2)(l 
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where is the energy of the electron/hole states referred 
to the Fermi level, E k = + A 2 is the quasiparticle 
energy, and A is the superconducting gap. 

The tunneling Hamiltonian, written in terms of the 
quasiparticle operators, is found by substituting Eq. (8) 
into Eq. (7). Eight of the resulting terms correspond to 
tunneling of a quasiparticle across the junction. For ex- 
ample, the tunneling from the state L0 to R0 is given 

by 



H T = t L R{u L u R e lip - v L v R e lip )~i L ^ m 
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where we have used thR = t* LR . As seen from the charge 
displacement operators, the two terms corresponds to 
charge tunneling in opposite directions. 

The tunneling of a quasiparticle produces both an 
electron-like and hole-like transfer of charge to the en- 
vironment. The superposition of charge transfer has not 
been considered previously for environmental theory, and 
can be incorporated readily into the formalism with only 
a change in the environmental displacement correlator 
(see Eq. (44) in Rcf. [1]) 
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where u = u L u R and v 



v L v R . The 4 correlators can 



be simplified by noting that for a harmonic oscillator the 
function is Gaussian, and therefore determined by its first 
and second moments [1]. Expanding terms to second 
order, one finds 
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where we have used (<p(t) — <p(0)) 
(<p(0) 2 ). One similarly finds 



and (<p(t) 2 ) 



(e 



¥(0>(0) 



(17) 



~ (2- 2tp(t)tp(0) - Mtf + v (0) 2 }) (18) 

= 2<l-[ V (t)+^(0)M0)) (19) 

~ 2e~ <[v(i)+v(0)lv(0)) . (20) 

Equation (13) can thus be written as 

e J(t) ={(u 2 + „ 2 ) e M*M°)> _ 2u ,, e -M*M0)>] 

x e -<*>(0)v(0)> . (2i) 



Note that the quantity (<p(i)(p(0)) is the noise correla- 
tion function for the (dimensionless) flux fluctuations. 
As shown in Eq. (2), it can be computed as the Fourier 
transform of the noise spectrum obtained from the fluc- 
tuation dissipation theorem. 

For the case Z r <C Rk, corresponding to large capaci- 
tance junctions where <p is small, the exponential in Eq. 
(21) can be expanded in a Taylor series 

e J « ~ (u - vf + (u + v) 2 (^(tM0)) . (22) 

Because the density of quasiparticles are small and they 
mostly thermalize to energies close to the gap, the oc- 
cupation probabilities are all approximately the same 
u l = v 1 = u r = v r = 1/2- In this limit, the proba- 
bility for a quasiparticle to tunnel without a change in 
energy is given by the (u — v) 2 ~ 0, showing a large 
suppression of tunneling due to the electron- and hole- 
like nature of the quasiparticle. However, the coherence 
factors when absorbing or emitting a photon to the en- 
vironment is proportional to the second term, which is 
near-unity u 2 + v 2 ~ 1. 

The coherence factors can be computed using Eqs. (9) 
and (10) 
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where in the second equation the terms £l£r have been 
averaged to zero for typical integrations over quasiparti- 
cle states. 

Wc first compute the quasiparticle tunneling rate (from 
left to right) corresponding to the first term in Eq. (22). 
This process is equivalent to a thermal current across the 
junction, as there is no charge transfer from tunneling. 
The average tunneling rate is 
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where Rt is the (normal state) resistance of the tunnel 
junction, and p = E/£ = Ej\jE 2 — A 2 is the normal- 
ized quasiparticle density of states. The factor of 2 for 
each p accounts for summing over both positive and neg- 
ative £. Note that without interaction to the environ- 
ment, the final energy of the quasiparticle excitation is 
unchanged from the initial energy. The integral, having 
only occupation factors, corresponds to simple tunneling 
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of the excitations through the junction, as appropriate 
for quasiparticles diffusing through the superconductor. 

The second term in Eq. (22) corresponds to a tunneling 
event that couples energy from the qubit to a quasipar- 
ticle. The rate for this process has been written down in 
the main article. 



QP DAMPING IN A COPLANAR RESONATOR 

In this section we compute the Q-factor (damping) for 
a coplanar resonator due to non-equilibrium quasiparti- 
cles in the superconductor. The discussion follows the 
Ph.D. thesis of J. Gao [2]. 

A bulk superconductor has complex conductivity a = 
u\ — i<j2 given by the Mattis-Bardeen theory 
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where parameters have been defined in the main text, 
and cr„ is the normal state conductivity We will assume 
that the quasiparticle density is small, giving o\ <C 02 
and little dissipation. 

The surface impedance of a superconductor depends 
on details of the superconductor and interface. In the 
thin film limit, the surface impedance is obtained via the 
geometrical factor of the film thickness d 
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The ratio of the real to imaginary impedance is 
Re(Z s 



Im(Z s ) 
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where here, for the thin film limit, 7=1. For the 
case of thick films in the local (dirty) limit, the surface 
impedance is given by standard electromagnetic theory 
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For o 1 <C (T2, the ratio of the impedances is the same 
form as Eq. (32), but with 7 = 1/2 in the local (dirty) 
limit and 1/3 in the extreme anomalous (clean) limit [3]. 

The quality factor Q may be computed taking into ac- 
count the field distribution in the coplanar transmission 
line. For a resonance frequency cj, standard microwave 
theory gives l/Q — TZ/luC, the ratio of the dissipative 
to total dispersive impedance. Here, £ = £ m + £ s and 
1Z = (g/s)Re(Z s ) are the total inductance and resistance 
per length of the line, £ m is the geometric inductance 
coming from the magnetic field, £ s = (g/s)lm(Z s /u) is 



the kinetic inductance, s is the width of the center line, 
and g ~ 1.0 — 1.2 is a geometric factor accounting for the 
coplanar geometry [4]. Defining the fractional contribu- 
tion of the kinetic inductance as a = £ S /C, one finds 
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For typical transmission lines the effect of the kinetic 
inductance is small, corresponding to a <C 1. 

In the thin- film limit, the fractional kinetic inductance 
using Eq. (31) is a ~ g/(ds<j2)/£ m , which gives 



1 



7r 2 V A / D(E F )A 



(35) 



where lZ n = \/dsa n is the normal state resistance per 
unit length of the coplanar line. The first term is the 
expected damping factor for a normal-state transmission 
line, while the remaining factors represents the suppres- 
sion of dissipation due to the small number of quasipar- 
ticles in the superconducting state. 

In the thick-film limit, the inductances may be param- 
eterized by lm(Z s /u) = fioX and C m = Hog m - In the lo- 
cal (dirty) limit, the penetration depth may be obtained 
from Eq. (33) 
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with a magnitude A <~ 50 nm for Al [5] . The magnetic in- 
ductance is described by a geometrical factor g m ~ 0.31 
that is only logarithmicly dependent on design parame- 
ters [4] . Using Eq. (34) , the quality factor is given by 
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Note that the Q scales as a n , which is a smaller power 
law than for the thin-film limit. 

We have not yet accounted for the change <5er 2 from 
n qp in Eqs. (29) and (30), which slightly decreases the ki- 
netic inductance and increases the resonance frequency. 
The fractional change in frequency is computed as 
Slu/uj = —S£ s /2£ = —a6£ s /2£ s = —a^5<j2/2a2 = 
a"f 8a\ /<72 = l/Q, so that the fractional frequency change 
is the same magnitude as l/Q. 
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